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C\j ! Abstract 

The fermion determinant in an instanton background for a quark field of arbitrary mass is 

""xl" ! 

studied using the Schwinger proper-time representation with WKB scattering phase shifts for the 

i-Q ■ 

relevant partial-wave differential operators. Previously, results have been obtained only for the 
extreme small and large quark mass limits, not for intermediate interpolating mass values. We 
show that consistent renormalization and large-mass asymptotics requires up to third-order in the 

• i— i , 

^ . WKB approximation. This procedure leads to an almost analytic answer, requiring only modest 

h : 

numerical approximation, and yields excellent agreement with the well-known extreme small and 
large mass limits. We estimate that it differs from the exact answer by no more than 6% for 
generic mass values. In the philosophy of the derivative expansion the same amplitude is then 
studied using a Heisenberg-Euler-type effective action, and the leading order approximation gives 
a surprisingly accurate answer for all masses. 
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I. INTRODUCTION 



The one-loop effective action plays a central role in quantum field theory. For a specific 
background field it represents quantum effects of direct physical relevance (e.g., the induced 
vacuum energy in the given background), while for general background fields it corresponds 
to the proper-vertex generating functional of the theory. The ultraviolet divergence or renor- 
malization problem for the one-loop effective action is well-understood, but it is notoriously 
difficult to evaluate explicitly its full finite part in the presence of a given nontrivial back- 
ground. For nonabelian theories in three spatial dimensions, there are essentially only two 
types of exact calculations for this full amplitude - the Heisenberg-Euler-type nonlinear 
action describing vacuum polarization and pair creation effects in the background of (co- 
variant*) constant gange held strengths flSflQ, and the QCD one-,„„p effective action 



for massless quarks in a self-dual Yang-Mills background field (describing instantons) [5|, | 

To study instanton-related physics it is of fundamental importance to determine the one- 
loop tunnelling amplitude given by the Euclidean one-loop effective action in the background 
of a single (anti-)instanton 7]. 't Hooft |5j succeeded in calculating the corresponding contri- 
bution by massless scalar or quark fields exactly in the 1970s, but such an exact calculation 
is no longer possible if the fields have non-zero mass. For various phenomenological applica- 
tions [8] (and also for the extrapolation of lattice results |9|, obtained at unphysically large 
quark masses, to lower physical masses), it is important to have more definite information 
about the contributions due to quarks of not-so-small mass. In this work, we will describe 
our approach to this problem and present explicit numerical results. 

Previously, two of us (with Kwon) studied the corresponding one-loop effective action 
with arbitrary mass, using a smooth interpolation between the results obtained in the large 
and small mass limits. The expression in the large mass limit is naturally obtained from the 
Schwinger-DeWitt expansion 12, (or the heat-kernel expansion) within the proper- 
time representation of the effective action, while the result in the small mass limit follows 
from the works of 't Hooft Carlitz and Creamer and Kwon et al. Q|. In this work, 
we use a systematic approximation for any quark mass - the effective action for any mass 
is calculated approximately without invoking an ad hoc interpolation procedure. The basic 
strategy is as follows. We use the proper-time representation of the effective action, which 
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requires the explicit functional form of 

F(s) = [ <fx tr(x|e~ s( - D2) - e~ s{ - d2) \x) (1.1) 



The instanton background enters through the covariant derivative D., in D 2 = D..D,,, and 

n n n 

this function F(s) can be expressed in terms of scattering phase shifts [5|, |15|, |16[ for the 
partial- wave quadratic differential operators related to D 2 . We evaluate these phase shifts in 
the quantum- mechanical WKB approximation, extended to the third-order correction terms 
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to ensure the correct small- s behavior for F(s). This leads to a numerical expression 
for F(s), which in turn yields the effective action for arbitrary mass value. The resulting 
mass dependence is fully consistent with the conjecture made by Kwon et al jl(3 ]. on the 
basis of the explicitly known results at the opposite ends. 

Our partial-wave WKB phase-shift method should provide a practical approximation 
scheme for the one-loop effective action in a broad class of background fields. In this ap- 
proach, F(s) generally takes a local form in the potentials and their derivatives - see (|3.11j) 
below. In the light of this observation, we conclude this paper by comparing to another 
method, the derivative expansion, which produces such local expressions. Surprisingly, the 
leading order of the derivative expansion provides a remarkably good approximation for 
general values of the quark mass. 



II. EFFECTIVE ACTION, PROPER-TIME REPRESENTATION, AND PHASE 
SHIFTS 

Due to a hidden supersymmetry of the system with a quark in a background instanton 
field, the one-loop effective action of a Dirac spinor field of mass m (and isospin |), T F (A; m), 
can be related to the corresponding scalar effective action (for a complex scalar of mass m 
and isospin |) by 0, 0, Q 

T F (A; m) = -~ In (™^\ - 2 V s (A; m), (2.1) 

The first contribution corresponds to the existence of a zero eigenvalue in the spectrum of 
the Dirac operator for a single instanton background. This relationship (which is special 
to a self-dual background) has the important consequence that it is sufficient to consider 
the scalar effective action V s [A; m) to learn also about the corresponding fermion effective 
T F (A; m), for any mass value m. 
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We consider an SU(2) single instanton background 
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A^x) = A^x)- 
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V^uaT Xy ^ _ 1,2,3,4; r = Jx^x^) 
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The regularized one-loop scalar effective action has the proper-time representation 

r°° ds 



(2.2) 
(2.3) 
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(2.4) 



where D 2 = D^D^ with = <9 M — iA^x). From this one obtains the renormalized effective 
action, in the minimal subtraction scheme, as 



V s {A- m) 



lim 

A^oo 



= lim 

A^oo 



T S A (A; m ) - 1 JL In (^j / d*x tr(F^) 



n(A; m) - 1 In f ^ 



(2.5) 
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Moreover, by dimensional considerations, we may introduce the modified scalar effective 
action V s (mp), which is a function of rap only, defined by 

1 



V s (A; m) = -ln(/ip) + f 5 (mp) 
6 



(2.6) 



and concentrate on studying the mp dependence of V s (mp). Then there is no loss of gener- 
ality in our setting the instanton scale p = 1 henceforth. 

The small- s behavior of F(s), as given by the Schwinger-DeWitt expansion, reads [lfll | 
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(2.7) 



Insering this result into (|2.4jl gives rise to the following large-mass expansion of V s (m): 
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(2.8) 



Note that the Schwinger-DeWitt expansion is a small s expansion for F(s), which naturally 
leads to a large m expansion for V s (m). To obtain an expression for T (m) for general values 
of m we need a more general expression for F(s). In the small mass limit, on the other hand, 
a completely independent calculation for f s (m) has been given (5, 10, based on the fact 
that the massless propagators in an instanton background are known in closed-form: 



m 







V s (m) = a Q^J + -{him + 7 - ln2)m 2 + 0(m 4 ), 



(2.9) 



where a(h) — 0.145873, and 7 ~ 0.5772 ... is Euler's constant. From this small m behavior 



it is possible to infer that F(s) ~ — j- as s —>■ 00. However, this information does not 
determine the magnitude of V s (m) even for small m, since the integral in (|2.4jl is dominated 
by contributions from non-asymptotic s-values. 

In this paper we make use of the fact that the function F(s) may be expressed in terms of 
appropriate scattering phase shifts. Note that the differential operator — D 2 in the instanton 
background (J2.2)) (with p — 1) can be cast in the form \^ 

d 2 3 d 4 4 , 4 



D 



dr 2 



3d 4 ?2 4 



J 2 



(r 2 + 1) : 



(2.10) 
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with L a — ^r] lJiUa x f j i d u (satisfying angular-momentum commutation relations) and eigenvalues 
L 2 l(l + 1), I = 0, ^, 1, -, ■ • ■ ; J 2 = (L + f) 2 with eigenvalues J 2 = + , j 

I I ± — I . Without any background, we have the differential operator 
2 



where T 2 = T a T a , and eigenvalue T 2 = t(t + 1)- appropriate to isospin t = -; L 2 = I/ a L a 



a 2 



9 2 



Sr 2 



-- + -L 2 , 

7" C^f i y >< ^' 



(2.11) 



which corresponds to the t=0 case of the expression (|2.1()jl . We may then consider the 
quantum mechanical scattering problem with the Hamiltonian 7i = —D 2 , viz., 



d 2 3d_ 41(1 + 1) 4{j + 1) 



dr 2 r dr 



r 2 + 1 



(r 2 + iy 



with the corresponding free Schrodinger equation given by 
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^0 = fc 2 *o. 



(2.12) 



(2-13) 



As r — > 0, we assume that \l/o ~ (const.) r 2 '. Also, to make the spectrum discrete, it 
is convenient to put the system in a large spherical box of radius R, demanding a suitable 
boundary condition at r = R (e.g., the Dirichlet condition i)j(r = R) = 0). Then the 
solutions of ()2.12j) and 1)2.13)1 have the asymptotic large-r behaviors 



#0n(r) ~ 2Cr~ 3/2 cos[k (n)(r + a)}, 
^„{r) ~ 2CV~ 3/2 cos[A;o(n)(r + a) + rj(k(n))] 



(2.14) 
(2.15) 



where rj(k(n)) denotes the related scattering phase shift, and the now discrete momenta 
k (n), k(n) (n: nonnegative integers) satisfy the conditions 1(| 



k (n + i) _ k (n) = ^ + O (± 



k (n + 1) - fc (n)), 



*b(n) = + + O . (2.16) 

This scattering mode description may be considered for every partial wave. If [k l,J (n)} 2 
and [&o(n)] 2 denote the energy eigenvalues introduced in association with ()2.12|) and (|2.13jh 
respectively, the function F(s) (see p.ljl ) may then be represented as 

F ( s )= E J2( 2l + 1 )( 2 J + 1 )J2{ e ^ [kl ' Hn)]2 -e- s[fc oHl 2 }, (2.17) 

including the degeneracy factor (21 + 1) (2j + 1). The factor (2j + 1) corresponds to different 
eigenvalues of J3, while the factor (21 + 1) corresponds to the eigenvalues of L 3 , the third 
component of the second set of conserved angular-momentum L a = — ^f] fiua x fl d u . 
The phase-shift relations (j2.1fij) imply that for large R we can write 



R S + °{W)\- (218) 

Based on this observation, it is possible to replace the sum J2 n i n (|^-17|1 by an integral: 

2 r°° , 

F(s) = -s £ dke- ks k(2l + l)(2j + l) Vu (k). (2.19) 

1=0,1,- ^ Jo 

Note that this expression is not only an infinite series but also contains an improper k- 
integral - hence it must be considered carefully. In the instanton background in particular, 
the /-sum and j-sum in (|2.19J) may not be considered in a completely independent way. This 
follows from the nature of the scattering problem as defined by (|2.12j) and (J2.13|) ; according 
to the forms for 7i and 7io, their small-r behaviors match for a given /-value, but it is the 
j-value that governs the large-r behavior of the effective potential in TC, and j does not 
appear in 7io- This apparent mis-match can be resolved simply by considering the phase 
shifts r] ll+ i(k) and 7] l+ ij(k) (with the same associated degeneracy factor) together as a 
package. With this understanding, the expression (J2.19|) can now be cast in the form 

F(s) = -s £ r^e- fe2s A;(2/ + l)(2/ + 2){r ?M+ i(fc) + ^ +i/ (A ; )}. (2.20) 

This form still requires careful treatment with regards to the /-sum and /^-integration. We 
will argue below that the correct (gauge- invariant) procedure is to have the various terms 
corresponding to the same 'energy' eigenvalue (i.e., the same k 2 - value) receive uniform con- 
sideration. As a convenient check on this procedure, we confirm that the predicted small-s 
behavior reproduces the form in (|2.7J) . 



III. SYSTEMATIC WKB PHASE-SHIFT ANALYSIS 



To find the exact form of F(s) with the help of ([2.20)1 . one must first have complete 
knowledge of the scattering phase shifts r)ij(k), and then carry out the needed infinite 
sum/integration in a carefully controlled way. This is not possible in general, and therefore 
one needs to develop a reliable approximation scheme to determine the function F(s). We 
provide such a scheme below, which relies on a systematic WKB approximation for the 
scattering phase shifts in question. In contrast to previous applications of the WKB 
method for a similar purpose (but in lower dimensions), the leading-order WKB result turns 
out to be insufficient even to ensure the correct value for F(s = 0). [Recall that from ([2.7)1 
we must have F(s = 0) = — ^ in order to construct the renormalized effective action as in 
([2.5)1 ]. The necessity to include higher-order WKB contributions arises from the inaccuracy 
of the /-sum in ([2.20)1 with respect to large-/ contributions. While the degree of accuracy of 
the WKB expression for the phase shifts is generally enhanced for higher partial waves, the 
large degeneracy factor (2/ + l)(2/ + 2) has the consequence that leading-order WKB method 
is insufficient. Fortunately, as we shall see below, this can be remedied in a systematic way 
by including higher-order WKB correction terms. 

Let us first consider the representation of F(s) in the leading WKB approximation. To 
that end, ([2.12)1 may be rewritten as an equation for = r 3//2 \I / (r): 
I 8> + m^+i + 40 - 00 + ! + !) _ 31 = 



I (^1T ^ >y* 2 iy* 2 j ^ 2 j ^ ^ 2 

For this one-dimensional Schrodinger-type equation, the leading WKB solution can be given 
immediately, including the usual Langer correction 19] to take into account the singular 
centrifugal term. Phase shifts in the leading WKB approximation, which can be extracted 
from this solution, read 



r${k) = r dr'Jk 2 - V l d{r') - [°° dr'yjk 2 - V£(r') (3.2) 



with 



(3.3) 



4(/ + 



1\2 



Vi(r) = ( 3 - 4 ) 



r 



In ([3.2)1 . r±(k), ro(k) denote the classical turning points determined by the conditions 
V l) i{ri) = k 2 and V^(r ) = k 2 , respectively. 



We now define 

X$\a) = £dke-*'kr$(k) (3.5) 

as this kind of integral is relevant in the construction of F(s) in ([2.20)1 . Using the expression 
()3.2|) for the phase shift and changing the order of integrations, this quantity may then be 
rewritten as 



dr 



dk e~ k2s kJk 2 - V l >Hr) 



fci(r) 



fco(r) 



dk e- k2s k^jk 2 - V l {r) 



(3.6) 



where ki(r) = yV l ' J (r) and ko(r) = y Vq(t). The k- integration in (jH.fi J) can be carried out 
explicitly to give 



Thus, the leading WKB expression for F(s) is: 



7T 



^^*'U (2,+i)(2(+2) f" 



sy M+ S(r) 



-«Vg(r) 



(3.8) 



Higher-order WKB correction can also be included. For the Schrodinger equation (|3.1jl . 
one can derive the 2nd-order and 3rd-order WKB phase shifts (i.e., r]\-{k) and r]ij (k)), 
incorpora ting the Langer correction in an appropriate manner, along the line discussed in 
Refs. jl?! lis! Leaving the somewhat involved details of this derivation elsewhere j^j], 



we shall here only report the results for (s) and X^ j (s) (which are related to the phase 
shift contributions of respective order by the integral relation of ()3.5|) ): 
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1 1 ~d 2 V l > j (r) 

s s — 

4r 2 12 dr 2 



\ Ar 2 12 dr 2 



(3.9) 



x (3) = _V7T_ I- 

w 4s 3/ 2 y 



5s 2 

32r 4 



s 3 d 2 V l ' j (r) s 3 d 4 V 



7s A fd 2 V l 'H 



A8r 2 dr 2 



288 dr 4 



1440 



dr 2 



5s 2 



s 3 d 2 V \r) 



s 3 d 4 V l (r) 7s 4 fd 2 V^(r) 



32r 4 48r 2 dr 2 288 dr 4 1440 \ dr 2 ^ 
Inserting these results into ()2.20|) we obtain the 3rd-order WKB expression for F(s) 

F(s) WKB = FW(s) + F®(s) + F®(s) 
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(3.10) 



-^-= [°° dr Y (2l + l)(2I + 2)fe-^ 1+ 'Mtf^(r)-e-<%'( 



2 



(3.11) 



where the H(r) functions are local functions of the potentials in (J3.3)) and ()3.4j) : 



4r 2 32r 4 i 12 48r 2 / dr 2 



s 3 rf 4yU M 7s 4 /^y/j/^X 2 

+ 1 > ( 3 - 12 ) 



288 dr 4 1440 \ dr 2 



•s- 



5s 2 /s 2 s 3 \ ffiH 



OV 7 -T 4r 2 ^ 32r 4 I 12 ^ 48r 2 I rfr 2 



s 3 dX(r) | 7, 4 / jX(r) \ 2 (313) 



288 <2r 4 1440 \ dr 2 

It is important that in evaluating (|3.1H1 the /-sum be done first, and then the re- 
integration. This rule regarding what operations should be done first follows [20] by adopting 
a definition of F(s) as given by the infinite 'energy' cutoff limit (see the remark at the end of 
Sec. [n}. In fact, only with this procedure is the correct small- s behavior for F(s) found (see 
below) - this can be taken as further (a posteriori) evidence for our rule. We have verified 
explicitly that carrying out the r-integration for individual partial wave contributions first 
and then summing leads to incorrect expressions. The difference for F(s) is ^, which leads 
to a spurious quadratic divergence in the s-integral. 

Based on the form ()3.11|) . one can determine the function -F(s)wkb numerically. The 
infinite /-sum in this formula can be handled using the Euler-Maclaurin summation formula 



2l| (some care is needed for small-r values), and the result is a rapidly convergent series. 
Numerical integration with the resulting function of r can then be performed with very high 
accuracy. Figure ^ shows a comparison of the WKB approximations for F(s) with one, two 
and three terms in the WKB expansion. We plot the function F(s)wkb (i-e., including up to 
3rd-order WKB terms), together with the corresponding plots for Fi(s) = F^(s) (i.e., the 
result based on the leading WKB phase-shift expressions only) and F 2 (s) = F^(s) + F' 2 '(s) 
(i.e., including up to 2nd-order WKB corrections). The value we found for Fi(s) at s = is 
—■h, while the correct value is — ^ (see (|2.7|) V Both 2nd and 3rd order WKB give the correct 
value for F(0). Thus, as remarked earlier, the leading WKB result alone is not sufficient here 
even for the renormalization discussion. With the 2nd order WKB expression, F 2 (s), there is 
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FIG. 1: Plot of F(s), showing the 1st, 2nd and 3rd order WKB approximations. 

no problem with renormalization, but, for its first derivative at s = 0, we found F^O) = — ^j, 
which does not agree with the small s result from the Schwinger-DeWitt expansion in (j2.7J) . 
On the other hand, going to 3rd order WKB, -F(s)wkb gives the correct value for this first 
derivative, with F'(0)wkb — With the small-s behavior satisfactorily taken care of, 
we expect good agreement with the large mass behavior of the effective action, as can be 
confirmed from Figure 121 For very large s (which corresponds to the m — > limit), the 3rd 
order WKB expression F(s)wkb approaches zero faster than the true F(s). Nevertheless, 
we show below that this 3rd order expression gives an excellent approximation even in the 
extreme massless limit. 

Our function -F(s)wkb can be used to determine the one-loop effective action for arbitrary 
mass value. We insert -F(s)wkb into (|2.4jl . integrate over the proper-time s (numerically), 
and renormalize according to (J2.5|) . Extracting V s (m), as defined in (|2.6|) . we obtain the 
plot shown in Fig. |21 Also shown in Fig. |21 are the curves based on the inverse mass 
expansion ()2.8|) and small mass expansion ()2.9|) . Evidently, our WKB-based plot corresponds 
to a smooth interpolation of the latter two curves, as conjectured in Ref.[lfJ earlier. The 
agreement at large m is excellent. For m — 0, our WKB-based prediction gives rise to 
the value 0.137827, which is about 6% off from 't Hooft's exact value, a(|) ~ 0.145873. 
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FIG. 2: Plot of T s (m), comparing the 3rd order WKB result with the exact extreme large and 
small mass limits. 

[If F 2 (s) were used instead of -F(s)wkb for this calculation, the predicted value for r s (0) 
would be 0.158084, which is about 10% off from the exact value]. The discrepancy from the 
exact result is expected to be largest for m = 0; hence, we estimate that our WKB-based 
prediction of V s (to), for arbitrary mass m, is good to 6% accuracy. 

IV. FIELD-THEORETIC DERIVATIVE EXPANSION APPROACH 

The field-theoretic derivative expansion provides a quick and extremely simple estimate of 
the one-loop effective action in an instanton background for any value of the quark mass m, 
and we show here that even the leading order term gives surprisingly good agreement. The 
philosophy of the derivative expansion is to compute the one-loop effective Lagrangian for a 
covariantly constant background field, which can be done exactly, and then perturb around 
this constant background solution. The leading order derivative expansion approximation 
for the effective action is obtained by first taking the (exact) expression for the effective 
Lagrangian in a covariantly constant background, substituting the space-time dependent 
background, and then integrating over space-time. 

For an instanton background, which is self-dual, we should base our derivative expansion 
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approximation on a covariantly constant and self-dual background: one may then set (in 
a suitable gauge) |22| = F^n a T a , where the abelian field strength is self-dual, 



and n a is a unit vector in color-space. Comparing with the large p limit of the su(2) 
instanton field strength (|2.3jl . we identify F^n a = —^rj^ va . The exact one-loop scalar 
effective Lagrangian for a covariantly constant self-dual field is well-known (see, e.g., Eq. 
(2.11) in |23j|). Substituting the instanton form we obtain the derivative expansion (DE) 
approximation (here we set p — 1, as before) 



/'scalar 
^DE 



(47T) 



00 ds 



(1+r 2 ) 2 



1 

3 



12 s 



„2\2 



(4.1) 



The leading derivative expansion for the effective action is then obtained by integrating the 
effective Lagrangian ()4.1|) over space-time. Note that (|4.1|) has been renormalized on-shell, 
so that p = m in (|2.6|1 . To study both the large and small mass limits it is useful to express 
(J4.1|) in a different, but equivalent, form using the identity: 



00 du 



M J 



-2ku 



U 



sinh it 



1 + 



/■°° dxx , / 



2 1 2 



(4.2) 



Then the large mass expansion of the effective action is straightforward: 
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m 4 (l + r 2 ) 4 
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(4.3) 



(2Z + 2) (2/ + l)(2/)(4/ + 3) m 4 ' 

where are the Bernoulli numbers. Note that the large mass expansion f)4.3|) of the lead- 
ing derivative expansion approximation begins with -\, rather than the true behavior, 
because for the covariantly constant self-dual field the corresponding Schwinger-DeWitt co- 
efficient vanishes by group theoretic traces. 

To study the small mass expansion, we perform the r integral in ()4.3|) to obtain an exact 
integral representation of the leading derivative expansion approximation 
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FIG. 3: Plot of T s (m), comparing the leading derivative expansion approximation with the 3rd 
order WKB result and with the exact extreme large and small mass limits. 

It is now simple to expand each term for small mass to obtain the leading behavior 

T S DE (A; m) ~ 1 m(m) + (A _ J_ m ( 4 8) - C'(-l)) + ^ m 2 ln(m) + . . . , m -> (4.5) 

The agreement with the leading small mass behavior (J2.9|) of the exact result is quite remark- 
able. The coefficient | of the ln(m) term agrees, as it must by virtue of the /3-function. The 
constant term (jr — ^ ln(48) — ('(— ~ 0.14301 is only 2% away from 't Hooft's value of 
~ 0.14587, and the coefficient of the m 2 ln(m) term is ^ ~ 0.433, compared to Carlitz 
and Creamer's result of 0.5. Figure El shows a comparison of the associated function T (m) 
with the WKB result of the previous section and the exact large and small mass limits. The 
agreement is surprisingly good for such a crude approximation j^j]. 



V. CONCLUSION 



In this letter we have presented a computation of the fermion determinant in an instanton 
background for all values of the quark mass. Using 3rd order WKB approximation for 
the scattering phase shifts we obtained a result which interpolates very well between the 
known extreme large and small mass results. We expect even higher order WKB terms 
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to improve the accuracy further. Also, the agreement in the small mass regime might be 
improved by determining the large s asymptotic behavior of the proper-time function F(s) 



generalizing the method of Barzinsky and Mukhanov |25j]. Finally, we showed that a very 
crude leading order derivative expansion approximation based on a covariantly constant 
self-dual background leads to surprisingly good agreement. 
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